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We present a general construction of initial data for Einstein's equations containing an arbitrary 
number of black holes, each of which is instantaneously in equilibrium. Each black hole is taken to 
be a marginally trapped surface and plays the role of the inner boundary of the Cauchy surface. The 
black hole is taken to be instantaneously isolated if its outgoing null rays are shear-free. Starting 
from the choice of a conformal metric and the freely specifiable part of the extrinsic curvature in 
the bulk, we give a prescription for choosing the shape of the inner boundaries and the boundary 
conditions that must be imposed there. We show rigorously that with these choices, the resulting 
non-linear elliptic system always admits solutions. 
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I. INTRODUCTION 

In this paper we consider the problem of specifying 
quasi-equilibrium multi-black hole initial data for the 
vacuum Einstein equations. This problem is important 
as a starting point for the numerical simulations of bi- 
nary black hole spacetimes. The initial data should be 
such that the black holes are in equilibrium with time in- 
dependent horizon geometries. Furthermore, the entire 
geometry of the spatial slice should, in an appropriate 
sense, be in quasi-equilibrium with minimal spurious ra- 
diation content. If both these criteria are satisfied, and 
if the black holes are in a roughly circular orbit around 
each other, we could reasonably expect the calculated 
gravitational waveforms to correspond to observations by 
gravitational wave detectors. 

The problem of finding such initial data has received 
a lot of attention in the numerical relativity literature 
and significant progress has been made in the last few 
years. The original numerical work addressing this issue 
was due to Cook in 2002 0| who, working in the so- 
called conformal thin sandwich (CTS) decomposition ^ 
of the initial data, proposed a set of conditions for solving 
the initial value problem subject to the quasi-equilibrium 
conditions. More recent developments in this direction 
can be found in 0, 0| . See Cook Q and Gourgoulhon et 
al. Q for an earlier numerical study of binary black 
hole initial data which is in quasi-equilibrium in the bulk. 
The case of quasi-equilibrium configurations inpresencc 
of matter is discussed in e.g. See also Q for an 

approach to this problem based on the post-Newtonian 
expansion. 

Independently of the numerical work, on the analytical 
side, a quasi-local approach to black ho le p hysics was de- 
veloped by Ashtekar e< aL (see e.g. [H El 111 El). This 
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work lead to the notion of an isolated horizon which mod- 
els a black hole in equilibrium in an otherwise dynamical 
spacetime. While isolated horizons are defined in the full 
four dimensional spacetime, it would seem natural that 
this framework should also have implications for the con- 
struction of quasi-equilibrium initial data. This issue was 
studied in some detail by Jaramillo et al. E3 ■ They used 
the isolated horizon formalism, also working in the CTS 
decomposition, to arrive at a set of boundary conditions 
for the constrained parameters of the initial data which 
is similar to that obtained by Cook, but with certain ad- 
ditional constraints on the boundary values of the other- 
wise free data. While the initial results are promising in 
both approaches, it is not yet settled if they will finally 
succeed. 

The aim of this work is to point out some potential 
mathematical difficulties with the approaches of 
and to suggest a resolution of these difficulties. Our aim 
is less ambitious than 0, 0, 0| in the sense that we only 
consider the issue of the appropriate boundary conditions 
at the horizon: i) we do not make any statement about 
the conditions that must be satisfied in the bulk to ensure 
that the entire data set is in quasi-equilibrium and ii) we 
do not discuss gauge conditions for the evolution equa- 
tions. Because we do not address i) and ii), we will work 
with the standard conformal transverse-traceless (CTT) 
decomposition of the initial data E3 ^^d not with the 
CTS as in P, 0, because the former simphfy the 
discussion. However we expect our main results to be 
relevant also in other decompositions including the CTS. 

The conceptual issues have to do with the boundary 
conditions for the momentum constraint. All the decom- 
positions involve an elliptic equation for a vector The 
references ^ 0, ^3 all use a Dirichlet condition for /3° 
which adapts the time evolution vector to the properties 
of the horizon. This requirement becomes intertwined 
with the initial data construction in the CTS decompo- 
sition since the latter can be interpreted to involve also 
the notion of time evolution. However, strictly speaking, 
the initial data construction is distinct from the choice 
of gauge for time evolution. More importantly, as al- 
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ready pointed out in a Dirichlet condition for 
entails potential problems for the very existence of solu- 
tions for the Hamiltonian constraint. We will show that 
a new type of geometric boundary condition for /J" can 
be used, which allows us to apply the theorems proved 
by Dain and Maxwell [Tfil in order to get a rigor- 
ous proof of existence of solutions for the resulting non- 
linear equations under appropriate assumptions. These 
assumptions, which involve restrictions on the shape of 
the inner boundaries and some inequalities the free data 
should satisfy, can be checked a priori. The advantage 
of this analysis is that the solution is guaranteed to exist 
in more general cases than the ones studied numerically 
so far; in particular, in non conformally flat cases. This 
is important for two reason: first, the conformal geome- 
try of a quasi-cquilibrium black hole initial data is still 
unknown (though it is clear that it cannot be confor- 
mally flat because the Kerr metric does not admit con- 
formally flat slices O). Second, one is usually interested 
not just in a single solution, but in families of solutions 
which depend smoothly on the relevant parameters of 
the problems, like separation distances, individual spins 
and linear momentum, etc. For example, one notion of 
quasi-equilibrium is defined as a variational problem for 
a family of solutions (cf. Q and also and references 
therein). Thus, we would like to ensure that the initial 
value equations admit solutions for the widest possible 
range of parameters and free data. 

Section^]sets up notation and states the mathematical 
problem we want to solve. Section IIIII gives the main 
result and compares it to earlier work, section llVl consists 
of a mathematical proof of the main result, and sectionlVl 
summarizes our results and suggests directions for future 
work. 



II. THE CONFORMAL METHOD AND 
NON-EXPANDING HORIZONS 

The problem we want to solve is to find initial data 
on a spatial slice M for the vacuum Einstein's equations. 
Thus, we want to find a 3-metric hat and a second fun- 
damental form 

kab ■■= -Kh^VcTd = -^Cr'hab , (2.1) 

such that the constraints are satisfied: 

baik"^ - kk'^) = 0, (2.2) 

n + k"^ -kabk'''' = 0. (2.3) 

Here t'^ is the unit timclikc normal to M, and Vq is the 
four dimensional covariant derivative. Our sign conven- 
tion for Kab corresponds to what is commonly used in the 
numerical relativity literature. Da is the derivative oper- 
ator compatible with hab, and TZ is its scalar curvature. 
We restrict ourselves to maximal slices, i.e. we always 
take K = 0. It is important to note that very little is 



known about solutions of the constraint equations (with 
or without inner boundaries) in the case when K is not 
nearly constant (see the recent review and references 
therein). 

The metric and the traceless part of the extrinsic cur- 
vature are conformally rescaled 

hab = i>^hab , k-' = V^'^A"^" . (2.4) 

As a rule, physical tensors on M will be denoted with a 
tilde to distinguish them from the corresponding confor- 
mally rescaled quantities (note the opposite convention 
with respect to Refs. [lllalli|l. In terms of the confor- 
mally rescaled quantities {hab, Kab), the constraint equa- 
tions become 

DaK'''' = 0, (2.5) 

Lhi' - -ii^afc^^"V"^ (2.6) 
o 

where L/j is the conformally invariant Laplacian operator: 
Lh ^ A - 7^/8; A := DaD" is the Laplacian, Da is 
the derivative operator compatible with hab, and TZ is 
its scalar curvature. To solve the momentum constraint, 
we decompose (the traceless) Kab according to York's 
prescription [l5l |: 

Kab = imab ~ Qab, (2.7) 

where /3° is a vector field on M, C is the conformal Killing 
operator defined as 

{CP)ab = ID^^aPb) - \habD,P'^, (2.8) 

and Qab is a freely specifiable symmetric and traceless 
tensor. The decomposition of Kab given by equation 12. 7() 
is known as the conformal transverse-tracelcss decompo- 
sition (CTT) but it is not the only possibility. In fact, the 
currently more commonly used decomposition in numer- 
ical relativity is the so called conformal thin-sandwich 
(CTS) decomposition also proposed originally by York 

Using the CTT decomposition, the momentum con- 
straint (|2.5|) becomes 

Al/3'' = J". (2.9) 

Here we have defined 
A^/3° Db{CI3Y^ = ^13" + ^D^DbfS^ + UIP^ (2.10) 

and J" := DbQ'^^- Thus, we have to solve the system of 
elliptic equations H2.6|l and (|2.9|) on a domain A/ C M"^ 
subject to certain boundary conditions which will occupy 
us for the rest of this paper. 

We wish to solve equations (|2.6() and (|2.9|l on a do- 
main M which is with an arbitrary finite number of 
compact sets excised from it. The inner boundary will 
be denoted by S, and in general, it is allowed to have 
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FIG. 1: The inner boundary S of the Cauchy surface M. In 
this figure the inner boundary consists of two disconnected 
spherical surfaces. 



many disconnected components. The excision surface S 
serves as an inner boundary for our problem and it rep- 
resents the surface of the black hole. This is understood 
to mean that each connected component of S is an outer 
marginally trapped surface, i.e. when it is viewed as a 
closed two surface in the full four-dimensional spacetime 
manifold, the congruence of future directed outgoing null 
geodesies orthogonal to it has zero expansion. This is de- 
picted in figure^which shows the Cauchy surface M with 
inner boundary S. The unit spacelike normal to S (with 
respect to the physical metric hat and pointing towards 
spatial infinity) is s°, the unit timclike normal to M is 
r", and 

r:=T" + S" (2.11) 

is a fiducial outgoing null normal to S. The ingoing null 
normal to S is 



The conformally rescaled normal s° is defined as 



-2 a 



(2.12) 



(2.13) 



The induced physical 2-metric on S is denoted by qab and 
it is conformally rescaled as 



qab (lab- 



(2.14) 



The projection operator onto S is = q'^. The second 
fundamental form of the inner boundary is defined by 



kab = qlq^DcSd- 



(2.15) 



We denote by k the trace and by the trace-free part 
of kab, that is 



k 1.-1° 

l^ab — 



b "I" "^kqab, k — q°'^kab 



(2.16) 



With this notation, wc can write down the expression 
for the expansion 0+ of 



e+ = r'^aib 



In terms of the conformally rescaled fields and using 
again = 0, this is equivalent to 



(2.18) 



Our first requirement on S is that the expansion 0+ 
must vanish. This is also known as the apparent hori- 
zon boundary condition in the literature. The problem 
of constructing initial data satis fying th is boundary con- 
dition has already been studied |1(t[ I17| . This boundary 
condition does not place any restrictions on the physical 
parameters of the black hole. We wish to impose the con- 
dition that the black hole is isolated in a certain specific 
sense. 

Our notion of S bein g is olated is based on the isolated 
horizon framework [lOl llll Il3. l2l| . This is a quasi- local 
framework to study black holes without relying on the 
global notion of an event horizon; see for a general 
review. It has found many applications in black hole me- 
chanics, mathematical relativit y, q uantum gravity, and 
also in numerical relativity [ijj Ull ■ An isolated horizon 
can be viewed as the world tube T of marginally trapped 
surfaces obtained by time evolution, in the case when T 
is null. In this case, it can be shown that the area of 
the black holes is constant in time and also the flux of 
gra vitational waves falling into the black hole vanishes 
|22l |2^ 0]. In the general case, T is expected to be 
spacelike, thereby representing a dynamical black hole. 
This is described by a dynamical horizon [23l [23 | . Trap- 
ping horizons [23. |25| can describe the null and spacelike 
cases in a unified framework. 

For our purposes, we do not need the details regarding 
isolated horizons. We only need to know that in vacuum, 
the condition for the world tube being null is equivalent 
to the vanishing of the shear cjab of pit |2^ . This is 
essentially a consequence of the Raychaudhuri equation 
whose expression, for a future directed geodesic -^^ with 
affine parameter A, is 
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Rc 



(2.19) 



KabS^S^ -K + k. 



(2.17) 



Here LOab is the twist of and Rab is the spacetime Ricci 
tensor; note that Rab^"^^^ = in vacuum. If 0+ is ini- 
tially zero, and taking the null rays to be surface forming 
{iJab = 0), then a non-zero shear would imply that 0+ 
will be non-zero at a later time. Thus the apparent hori- 
zon cannot evolve along ^°^ , and the black hole will not 
be isolated; see |23,|3| for further details and discussion. 
Conversely, if we take S to evolve strictly along f " so that 
T is null and 0+ is zero at all times, we get the condition 
that the shear vanishes identically, which in turn can be 
shown to imply that the rate of change of the area of 
S is zero, and S can be considered to be isolated. In 
the presence of matter, we must impose, say, the null en- 
ergy condition to get analogous results. The condition 
(Jab = has also been considered independently by Cook 
[l|. In the language of the isolated horizon framework, 
this condition would guarantee that S" is a cross-section 
of an infinitesimal non-expanding horizon [T]| . 
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The shear of £° is defined as the tracefree part of the 
projection of V(a^5) onto S: 

<yab = qlqt^iJd) - ^e+g,b . (2.20) 
In terms of the physical fields it is given by 

Oai, = -qlqtKcd + Iqatq^'Kcd + ~k°b , (2-21) 

and. in terms of the conformally rcscalcd fields, the shear 
is written as 

<yab = [-q'aqtKcd + Iqabq^'^K^a^ + i^^kl, . (2.22) 

Here kab is the extrinsic curvature of S embedded in Ad 
with respect to the conformal metric hab 

kab = qaqtDcSd, (2.23) 

and, as before, 

kab = C + ^kqab, k = q°-^kab. (2.24) 

Our task is to prescribe appropriate boundary condi- 
tions for (■0,/3") so that equations H2.6(l and 1)2.9(1 have a 
regular solution with everywhere positive, and = 0, 
CTah = at 5. 



III. THE MAIN RESULT 

A. The sign of i^abs^s' 

Let us start with the condition 0+ = which, from 
equation H2.18|l . can be written as 

As-'Dai' = -i^k - i^-^Kabs'^s'' . (3.1) 

We would like to use this as the boundary condition for 
solving equation H2.6|) for -0, together with the standard 
condition at infinity 

lim V = 1. (3.2) 

r — *oo 

Let us first discuss some general properties of Eqs. ((2.6|l . 
(Tni) and (E3. 

We begin with the physical meaning of the sign of 
Kabs"^s^. In a realistic collapse, the boundary is ex- 
pected to not only satisfy 8+ = 0, but also to be a future 
marginally trapped surface, that is 0_ < 0, where 0_ 
is the expansion of the ingoing null-normal n° defined in 
Eq. I|2.12|l . For example, a surface with 6+ = and 
0_ > is not expected to be present in a realistic grav- 
itational collapse of matter; these surfaces are located 
on the inner null boundary of the left quadrant of the 



Kruskal diagram (region IV of Fig. 6.9 in H^). The 

expansion 0_ is given by 

e_ = -k + kabs°-s^ - k (3.3) 
= i'-^ {-'is'^Dai' -^k + ^r'^Kabs'^s'') , (3.4) 

where the second line applies if the data is maximal {K — 
0). 

Under this assumption, we get that 8-|- -|- 0_ = 
2KabS°^s'^ 1 then for a future marginally trapped surface 
on a maximal slice we always have 

/^a6s"s' < 0. (3.5) 

Also, as mentioned below, it turns out that we also need 
to control the sign of KabS°'s^ in order to provide nec- 
essary conditions for the existence of solutions of Eqs. 
(E21), EU and 

On the other hand, if we were to impose purely Dirich- 
let boundary conditions on , then we could not guaran- 
tee that Kabs'^s'' has a definite sign at the boundary. This 
is simply because equation H2.7|l implies that KabS°'s^ in- 
volves radial derivatives of /3° which cannot be controlled 
by Dirichlet conditions: 

KabS'^s'' = (/:/?)ahS°s'' " Qabs'^s'' . (3.6) 

This conclusion is also valid for the conformal thin- 
sandwich decomposition. This is because also in the GTS 
decomposition, the extrinsic curvature is written in terms 
of the derivative of Z?*^ and KabS°'s^ will again involve ra- 
dial derivatives of Z?*^. 

We would like to prescribe [CP)abS°'s^ as free data on 
S and, at the same time, be able to enforce Uab = on 
S. The main result of this paper is that this is indeed 
possible. 



B. Existence of solutions to the Hamiltonian 
constraint 

Let us now review the condition telling us about the 
existence of solutions to equation (|2.6|) with boundary 
conditions 1)3. l|l and H3.2|l . The conformal factor i/' ap- 
pears in the denominator of both equations H3.1|l and 
(|2.6|) . hence these equations are singular if the confor- 
mal factor vanishes, which is in agreement with the fact 
that the physical metric hab does not make sense if the 
conformal factor is zero at some point. For general con- 
formal metrics and boundaries, these equations will have 
no solutions. To illustrate this, let us consider the time 
symmetric case [Kab = 0). For the case with no in- 
ner boundaries, necessary and sufficient conditions on the 
conformal metric hab for the solvability of Eq. (|2.6|) and 
(|3.2|l has been studied in [23| and 0|. This condition 
can be written in term of the following quantity known 
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as the Yamabe invariant 



inf 



Mile 



(3.7) 



where C^{M) denotes functions with compact support 
in M, WfWle = (/m Ifl'^y^^, and ^ means that ip 
cannot be identicaUy zero everywhere. The number Ah 
is conformally invariant in the following sense: let 9 be 
any smooth and positive function, and let hab = dhab', 
then = Xh- The Yamabe invariant has a long history. 
It was discovered by Yamabe for compact manifolds and 
later studied by many authors; see, for example, the re- 
view m. 

A solution of (|2.6|l - (|3.2() (with Kat = 0) exists if and 
only if Xfi > 0. Every metric with TZ > (the flat metric, 
for example) satisfies this condition. It is obvious from 
(|3.7() that 7^ > implies Ah > 0, however, to prove that 
in fact A/i > is non trivial. However, there exist metrics 
which do not satisfy Ah > (see for an explicit, ax- 
isymmetric example). For these conformal metrics there 
are no solutions of (|2.6() - (|3.2() . 

If we include inner boundaries, there exists a general- 
ization of this condition (cf. ^^) in terms of the follow- 
ing conformal invariant which is a generalization of the 
Yamabe invariant: 



Ah,s = inf 



(3.8) 

This invariant has been studied also for the compact case 
in [23. We have the following result proved in 0|. A 
solution of Eqs. ^i^, and ||22J), with R'ab = 0, 

exists if and only if Ah.s > 0. Note that now Ah,s depends 
on the choice of the boundary. In particular, note that 
the boundary term in H3.8|) has a minus sign.^ It can be 
proved that any metric with TZ > and boundary with 
A: < have Ah,s > |0. As before, in this case it follows 
directly from the definition that Ah,s > 0, but to prove 
that it is strictly positive requires extra work. Note that 
this does not apply to the flat metric with spheres as 
boundaries, since for spheres of radius ro we have k = 
2/ro > 0. However, the flat metric with an arbitrary 
number of non-intersecting spheres as inner boundaries 



^ In this article we always assume the dimension of M is 3, however 
all of the follo wing discussion can be generalized to arbitrary 
dimensions, see lid . 

^ Unfortunately, there exists in the literature different conventions 
for the signs of K^b, k, and s". In this article we have used what 
seems to be the standard conventions in numerical relativity. The 
relation of our present convention with the ones in Il6l and llTl 
is the following. The second fundamental form K„i, is denoted 
by (Jab in lid , and let us denote by Kab the one in llTl . Then we 
have Kab = f^ab = ~Kab- The mean cur vatu re of the boun dary 
(in our notation k) is denoted by h in Il6l and by H in llTl . 
We h ave k = —2h = —H. The normals are denoted by i/" in 
both and ^3 (they use the same choice for it). We have 



satisfies Ah,s > 0, provided they are not too close to each 
other. This can be seen as follows. Assume that we have 
only one sphere of radius ro centered at the origin. Take 
the conformal factor 6 = 1 + ro/r, and define the rescaled 
metric hab = ^'^^ab (where Sab is the flat metric). This 
is, of course, the Schwarzschild initial data. This metric 
satisfies R = and the boundary r = rg is a minimal 
surface with respect of the metric hab, i-e, fc = 0. From 
the previous discussion it then follows that A,- 5 > 0, and 
since it is conformally invariant we have Xs^s — Xf^ g > 0. 
In order to generalize for more spheres we will use the 
well known Misner initial data: Misner shows in Lemma 
3 of I^Ol that there exists a conformal factor 9 such that 
hab = O'^Sab, R = 0, with the boundaries of the spheres 
being minimal surfaces with respect to the metric hab, i-e. 



k = 0. Then we have also in this case A^^c 
For two spheres, the condition that the spheres are not 
too close used in this Lemma, just means that they do 
not touch each other. 

It is easy to construct more general examples of metrics 
which satisfy Xh,s > 0- As pointed out in [l^, if we take 
any metric on R"^ which satisfies TZ > and we excise 
appropriate small spheres on it, we get Ah,s > on M 
(see [i3 for details). 

For the non-time symmetric case, if we assume the 
maximal condition K = and no inner boundaries, then 
is clear that every maximal initial data satisfies Ah > 
because TZ = KabK"'^ > 0. This is also a sufficient con- 
dition (with appropriate fall-off conditions for Kab at in- 
finity) for the conformal metric to ensure the existence of 
solutions to the non linear equation (|2.6|) with boundary 
conditions H3.2|l (see [3j| and references therein). 

If we include inner boundaries, still requiring the con- 
dition K = 0, it is not clear if the condition Xhs > 
is necessary. However, it has been proved in 0, 0| 
that it is a sufficient condition together with additional 
assumptions on k and Kabs'^s''. Moreover in these ref- 
erences the condition Ah,s > plays an essential role. 
Even if the condition Ah,s > turns out not to be neces- 
sary in generic situations, it is very likely that it will still 
be relevant for the black hole problem. The reason being 
that for these data, time symmetry arises as a limit when 
the linear and angular momentum of the black holes are 
zero, and in this discussed above, Ah,s > is a 

necessary and sufficient condition. 

The assumptions on k and Kabs"'s'' are always made in 
some representative metric in the class Ah,s > 0- That 
is, these assumptions are not conformally invariant. The 
two main examples are the following. If we assume that 
the conformal metric satisfies TZ> and fc < (this auto- 
matically implies Ah,s > 0) and we assume Kabs'^s^ > 
at the boundary S, then there is always a solution of Eqs. 

(|TT|) and lIX^ (with K = and appropriate fall-off 
conditions for Kab)- This case can be obtained from [l7j 
making a time reversion (that is, t°- — i", 8+ — 
0_ — > — 0+, Kab ^ ^Kab)- This example seems to be 
physically relevant only if KabS°'s^ = at the boundary 
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since, if this is not the case, then 0_ > 0. The second 
example was studied in [T^ . The conformal metric is as- 
sumed to satisfy Xh,s > and, in addition, TZ = 0, k > 
(the flat metric with spheres as boundaries satisfies these 
conditions). It is assumed also that —k < KabS°'s^ < 0. 
In this case we get future trapped surfaces. 

It is important to note that for an arbitrary metric 
which satisfies Xh.s > 0, it is possible to calculate (solv- 
ing a linear equation) a conformal factor 6 such that the 
rescalcd metric hab ~ d'^hab satisfies R = and fc > or 
(with, of course, a different 9) R> and k <0 (cf. [lfil|V 

We have seen that in both cases the existence theo- 
rems assume a definite sign for Kabs'^s^ at the boundary. 
This is essential in order to control the positivity of the 
conformal factor using the boundary condition and 
the maximum principle. In order to be able to control 
KabS°'s^ at the boundary, it is necessary to use a bound- 
ary condition for the momentum constraint H2.5II that 
allows us to prescribe Kabs"s^ as free data. In (IQ, 0| 
this has been achieved using the following Neumann type 
boundary condition for /3° on S: 

{Cl3)abs''^Q. (3.9) 

However, we will see in the next section that with this 
boundary condition, we cannot enforce Gab = at the 
boundary. In order to do this a new kind of boundary 
condition will be needed. 



C. Boundary conditions for the momentum 
constraint 

Let us now focus on the condition (Tab ~ with <7ab 
given by equation H2.22|l . First note that the two terms in 
(|2.22(l involve different powers of the conformal factor ip. 
Since ip is the unknown solution of the nonlinear equation 
(|2.()|) . it will be very difficult to get aab = without 
requiring that both these terms vanish independently. In 
particular, this requirement will imply that the shear of 
the ingoing null geodesies along n° is also zero. Then, the 
second term of (|2.22f) suggests that we impose = 0. 
This is a condition on the shape of the boundary S, and 
is also known as the umbilical condition. For example 
a sphere in the flat metric satisfies this. Note that this 
condition is conformally invariant, that is fcj]^ = <S=> 

Having made the above choice for kab, we are then 
left with the first term of H2.22|l . Before discussing the 
general case, let us mention some important examples of 
initial data which not only satisfy <7ab = at S', but also 
have kab — 0. 

The first example is provided by the Misner solution, 
already discussed above. The fact that these data satisfy 
Uab = at S' can be directly verified from H2.21|l : due 
to the reflection isometry of these data we automatically 
have kab = |32l |. and Kab = because the data is time 
symmetric. 



The second important example is the Kerr initial data 
for the Boyer-Lindquist slicing. Since Kerr is stationary, 
it is clear that Oab = 0. Moreover, these data also have 
an isometry which leaves the horizon fixed, then we also 
have kab = in this case. Note that although the Kerr 
data are not time symmetric, the terms with Kab in (|2.21|) 
still vanish. 

Finally, some of the Bowcn-York data [s^ also satisfy 
these conditions. These are the data for one black hole 
with spin (defined by Eq. (10) in js^) and for the so 
called "negative inversion" single black hole case with 
linear momentum (defined by K~j in Eq. 9 of H^)- For 
these cases, we have an isometry which leaves the horizon 
fixed, so that kab ~ 0. Moreover, one can explicitly check 
that the first term in (|2.22|l vanishes for these conformal 
second fundamental forms. 

Let us now return to the general case. Wc make a 
decomposition of f3 into its normal and tangential parts 
with respect to S 

P'' = bs" + (3^^, (3.10) 

where /3f|sa = 0, 6 = f3''sa. We insert lETj) in fTT^ . we 
use the decomposition ()3.10() and after some computa- 
tions we find 

aab = (^' - 2^.-26)C 

+ (qMcd - \qabq''''Qcd - (//3||)afc) . (3.11) 

Here I is the conformal Killing operator on the surface S*, 
that is 

{lm)ab = d(a'mb) - qabdcm", (3.12) 

where d is the covariant derivative with respect to qab, 
is any tangential vector to S {s°'ma = 0) and = 
qabm''. 

The important point is that Eq. H3.11|l . if we as- 
sume fcj^jj = 0, only contains tangential derivatives of 
/Sll"^. These derivatives can be controlled using Dirichlet 
boundary conditions for Thus, if we choose (//3|| )ab 
and Qab appropriately, we can ensure the vanishing of 
the shear at S. It turns out that this is in fact a valid 
boundary condition for P"". 

The complete set of boundary conditions on (3°- for 
solving equation (|2.9|l is 

s'^s'imab = /, (3.13) 
= V'" , (3.14) 

where / and (y9° arc to be prescribed a priori, subject to 
certain restrictions mentioned below. 



^ Incidentally, this is also the case in 0, 1141 , where it is achieved 
by canceling the factor multiplying fcj^j^ by means of a Dirichlct 
boundary condition on b. The price is the loss of the control on 
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With these choices, the problem is solved as follows. 
Chose any and let Qab be such that 

{lv)ab " QaltQcd + \qabq'"'Qcd = 0, (3.15) 

(alternatively, one could prescribe Qab and attempt to 
solve for (^°). Solve equation 12. 9|) for using the 
boundary conditions l|3.13|l - 13.14|l . Then (provided 
^ah = 0) we will have aab = 0. Furthermore, KabS°'s^ = 
/ — Qabs'^s'^ is a free data which can be chosen to be 
non-positive by an appropriate choice of the quantities 
Qabs"'s'' and / which have, up to now, not been con- 
strained at all. The proof that it is always possible to 
solve equation (|2.9|1 with boundary conditions (|3.13|) - 
(|3.14|) is given in section Hvl 

Having solved the momentum constraint and having 
retained Kabs"'s^ as a free data, we can now solve the 
Hamiltonian constraint H2.6|l for the conformal factor 
subject to the boundary conditions H3.1|l and H3.2|) . The 
possible choices of this free data have been given ear- 
lier in section UlI CI Recall that to solve the momentum 
constraint, we only have fixed the tracefree part /cj^^; the 
trace k is still free. The physically most relevant case 
is when we want to obtain future trapped surfaces. As 
discussed in section UlI CI the two possibilities are 

i) 7^ > 0, /c < and Kabs'^s^ = 0. Note that aU 
the examples discussed above satisfy KabS°'s'' = 0. 
This condition implies 9+ = 0_ = on S. 

ii) The more general choice TZ — 0, k > 0, and KabS^'s^ 
such that 

-k < Kabs'^s'' <0. (3.16) 

In ii) the case Kabs"'s^ = is also included. However, in 
practice, for non conformally flat metrics, it is perhaps 
more convenient to work with i) because in this case the 
conditions 7^ > 0, /c < involve only inequalities which 
are easier to achieve (for example, by a small perturba- 
tion of a given metric) than the condition TZ = Q. 

Just like the Dirichlet and Neumann conditions, the 
boundary conditions (|3.13|) - H3.14|) are natural for the op- 
erator Ai. This is easily seen by considering the Green 
formula for the non flat operator A^,: 

/ {cpT^au = - / /3"AiC„ + i imabs'^e ■ 

Jn Jn Jan 

(3.17) 

Taking to be a solution of the homogeneous problem: 

AlC^O, (3.18) 
with boundary conditions 

s^'s'iCOab = 0, (3.19) 
?M = 0, (3.20) 



and /3° = the boundary term in equation H3.17|l can 
be seen to vanish. Thus we see immediately that 
must satisfy {C(,)ab = 0. Therefore, as for the Dirichlet 
and Neumann (i.e. equation 13.9|l ) cases, the kernel of 
the homogeneous problem with our new boundary con- 
ditions consists only of conformal Killing vectors. The 
only difference is that in the Dirichlet case the kernel 
is empty because there are no conformal Killing vectors 
which vanish at the boundary. In the Newmann case 
every conformal Killing vector is in the kernel because 
{C£,)abs'' vanishes identically. For our present conditions, 
the only conformal Killing vectors that are in the kernel 
are the ones which are normal to the boundary. This 
suggests that these boundary conditions are well posed, 
and we shall show in section Hvl that this is indeed so. 

Finally, we want to point out that the new boundary 
conditions (j^T^ - (jXT^ can be used in both references 
[TgI ITtI in replacement of H3.9|l , all the results presented 
there will hold without any essential change. 

D. The 2+1 decomposition 

While the above description is a complete specification 
of the problem, in order to get a better feeling of the 
nature of this elliptic system, we decompose the elliptic 
system into its radial and tangential parts, and show that 
the above boundary conditions translate into a Robin- 
type condition for the function b appearing in equation 
(|3.1U|) . Using the 2-1-1 decomposition I3.10|l of the shift 
to express the boundary conditions (|3.13() and (|3.14|l on 
the compact excised surface S, we find 

2s'' Dab -kb=^f + daV" - 2ip''da lu N , (3.21) 

= (^V (3.22) 

Here S is characterized as the inverse image of a constant 
value To by the height function r, and N is the normal- 
izing factor of such that Sa = NDar, and s°'Sa = 1. 

This decomposition suggests a general 24-1 splitting of 
A'/, provided by the slicing defined by r, as a procedure to 
solve the elliptical equation (|2.9|) for the shift. Extending 
(|TTU|) to the whole of M, we find for /3f|: 

Al/3|| = J"-Al(5s") = Ji'f . (3.23) 

We can enforce the radial part of J|° to vanish, if we 
impose the following elliptic equation on b 

Jl\Sa AL(5.s").Sa = rSa = J 1_ . (3.24) 

Expanding the differential operator we find 

J_L = DaD'^b + ^s^'s'^DaD^b (3.25) 

+ -Dabiks" - rf" In TV) 
3 

+ b (UacS'^s'' + ^D^'k - kack''' - {da lu N){d'' \llN)\ 
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Note that the principal part of the operator defined by 
(|3.25|) is not the Laplacian, but it is nevertheless elliptic. 
In this scheme, we calculate the scalar source from 
the original source J° by contracting with the normal 
vector s". Then we solve Eq. (|3.25() for b by imposing 
boundary condition (|3.21() . With the resulting b we cal- 
culate the source J|| and solve Eq. (|3.23|) for f3^^ by using 
the Dirichlet boundary condition (|3.22|) . 

We conclude with some remarks. First, there is a re- 
markable choice for namely to prescribe it as a con- 
formal Killing vector on S (this is condition (49) in Q 
or constraint (57) in ^i^sd is {l(p)ab = 0. This is a 

particular case of condition (|3.15|) . However, this choice 
seems to play no fundamental role in our analysis. 

Secondly, it is interesting to remark on a generalization 
of the above procedure. In 0, the free data is chosen in 
such a way that it is possible to control the size of both 
|8_| and |8+|. In other words, it is possible to control, 
in some sense, how trapped the boundary is. The same is 
possible here with the shear. This question is of interest 
if we want to control the amount of radiation falling into 
the black hole [H HI |23|. In equation (pTTs]) . if we 
choose Qab such that the right hand side is some Sah 
instead of zero, then the elliptic equations for and 
-0 can still be solved and the final solution will have an 
inner boundary with shear aab = — V'~^^ab- Since ^p > 1 
we get \aab\ < \'^ab\- It is also possible to get a lower 
bound for \aab\ using the upp er bound for tp obtained in 
the existence theorem in llTl. 



IV. ELLIPTIC BOUNDARY CONDITIONS FOR 
THE MOMENTUM CONSTRAINT 

It is well known that the operator defined by Eq. 
(|2.10() is elliptic. For a given elliptic operator, a set of 
boundary conditions are called elliptic if they satisfy the 
Lopatinski-Schapiro conditions, also known as the cover- 
ing conditions. For the definition of these conditions as 
well as the other concepts used in this section see, for ex- 
ample, the review |34j and references therein. These con- 
ditions are important because an elliptic operator with 
elliptic boundary conditions will always have solutions 
provided the sources and the boundary values satisfy a 
finite number of conditions. 

The operator A/, is (for a three dimensional manifold) 
of degree 3, and we therefore need to prescribe three 
equations as boundary conditions. Let O be a bounded 
domain in R"^ . The most important example of an elliptic 
boundary condition for A^ is the Dirichlet one 

13" ^ ip" on dn, (4.1) 

where if"" is an arbitrary vector on the boundary. This 
condition has been extensively used in numerical rela- 
tivity. The analog to the Neumann condition for A^ is 
given by 

iCf3)abs'' = fa on an, (4.2) 



where s° is the normal to dil. These conditions are el- 
liptic. 

In this section we want to prove that the following 
boundary conditions for A^ are also elliptic 

s''s''{C(3)ab = fi (4.3) 
Pami = ip2, (4.4) 
Pam'^ = ip-i, (4.5) 

where m° and m2 are tangential and linearly independent 
vectors at the boundary. Note that Eqs. H4.3|l - (|4.5|l are 
three linear conditions for /?" at the boundary. In order 
to write the Lopatinski-Schapiro conditions we need to 
define the principal part of the both the operator and 
the boundary conditions at an arbitrary point xq on the 
boundary. For the operator A^, the principal part is 
given by the standard definition. That is, it is given 
by the terms which contains only two derivatives. If we 
choose coordinates at xq such that hab{xo) = Sab, then 
the principal part of Aj^ is given by 

AO = AO/3'' + ^d''dbf3\ (4.6) 

where AO is the flat Laplacian and d denotes partial 
derivatives. 

For the boundary conditions we need to be careful in 
the definition of the principal part. If we choose only the 
terms which contain the highest order derivatives, then 
only Eq. H4.3|l will survive and this will lead to an ill 
posed problem. In order to take into account the fact that 
terms of mixed order appear in the boundary operator, 
we need to use the general definition of the principal part 
given by |35| . This definition involves integer weights 
that, with the notation of [s^l and for the operator A^, 
are given by = Sj^ = 1 [v = 1,2, 3). For the boundary 
condition, again in the notation of |3J|, we have tq = 0, 
ri ~ r2 = —1, where rp corresponds to Eq. (|4.3|) . ri 
to Eq. H4.4|l . and r2 to Eq. H4.5|l . With this choice the 
principal part will include all the equations (|4.3|I - H4.5|) . 
That is, the principal part at xq is given by 

s''s\C''(3)ab, f3am1, (3am^, (4.7) 

where 

{C°P)ab = 2d^al3b) - ^SabdcP', (4.8) 

and the tangential vectors m are evaluated at xq. 

For a given point xq at the boundary, we choose co- 
ordinates {xi, X2, X3) such that the normal is given by 
s = d/dx3, and {xi,X2) are coordinates on the tangen- 
tial plane at xq; see figure EJ Consider the homogeneous 
constant coefficient problem, on the half plane X3 < 
with boundary X3 = 

Alf]"" = (4.9) 

s''.s''iC°P)ab - (4.10) 
(3am.l = (4.11) 
Pam^ = 0. (4.12) 
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FIG. 2: The set-up of the Lopatinski-Schapho condition, xq 
is a point on the boundary dQ, {xi,X2) span tlie tangent plane 
to d^l at 2:0, and x^ is orthogonal to dO,; s" is the unit normal 
to (5f2. The shaded region is x^ < 0. 



The boundary conditions are said to satisfy the 
Lopatinski-Schapiro conditions if there are no nontrivial 
solution of H4.9|l - (|4.12|l of the form 



(4.13) 



where ^1 and ^2 are arbitrary real numbers, and v^{x3) 
tends to zero exponentially as 2:3 — > —00. To prove this, 
the key will be the following Green equation valid for any 

! Jn Jon 

(4.14) 

The proof is very similar to the one given in example 10 
of [SJ for the boundary conditions H4.1|) and (|4.2|) . Let 
us assume that there exists a solution of the form (|4.13() 
to equations (|4.9|I - H4.12|) in the half plane < 0. Let 
Li = 27r/^i and L2 ~ 27r/^2- Consider the following sub- 
domain contained in the half plane X3 < 0: the infinite 
cubic region X3 < 0, < iCi < Li, < 0:2 < ^2- For this 
sub-domain we use equation (|4.14|) for /3° = . We want 
to prove that, on this domain, the boundary integral in 
(|4.14|) vanishes for a solution of the form (|4.13|l . This 
is clear for the face — ^00 because the solution, by 
hypothesis, decay exponentially at infinity. Take the face 
X3 = 0. On this face the normal s° is also the normal to 
n. Eqs. imply that £°_ = as"" on this face 

for some function a. We use Eq. (|4.10() to conclude that 
the integrand in the boundary integral vanishes on this 
face. On the other faces, the integrand does not van- 
ish. However, because of the choice of Li and L2, we 



have that the integrand of opposite faces are identical. 
Then, the sum of the boundary integrals vanishes be- 
cause the normal is always outwards. We conclude that 
{C'^P)ab = 0. But there are no conformal Killing vec- 
tors in flat space which decay to zero at infinity. Hence 
the Lopatinski-Schapiro conditions arc satisfied, and the 
boundary conditions arc therefore elliptic. 

From the previous discussion, using standard results 
in elliptic theory, we deduce that a solution /3° of the 
boundary value problem 



j° on n 

f on dfl 
(pi, on dCl 
Lp2, on dfl. 



exists if and only if 



dn 



b.f= / 



(4.15) 
(4.16) 
(4.17) 
(4.18) 



(4.19) 



for all conformal Killing vectors of the metric hat 
which are normal to the boundary of $7, and where b 
is defined by = 6s". This can be shown by considering 
the Green equation H3.17|l from which wc deduce 



/ (rAi/3a - /^-^AiCa) 

Jn 



{{mabs'^e - mabs'-p') . (4.20) 

dn 

In this equation, set {C^)ab = (which implies A^,^" = 0) 
and = 6s" at the boundary, to immediately get Eq. 
(ICTll . 

If the metric admits no conformal Killing vectors, then 
there are no restrictions on J" . An example of a metric 
with conformal Killing vectors is the flat metric. If the 
boundary is a sphere centered at the origin, then = 
is a conformal Killing vector which is normal to the 
boundary (here we have assumed coordinates such that 

hab = 5ab)- 

In our case, J" = D^Q'^^ for some trace-free tensor 
Q°^^ . Then condition 14 . 1 91 can be written as 



an 



{bf - Qabs'C) = 0. 



(4.21) 



for the exterior region M discussed in this article, the 
boundary integral H4.21|l contains two terms, one is the 
inner boundary S and the other is an integral over the 
sphere at infinity. 

For simplicity we have considered in this section only 
bounded domains, since the new part here is given by the 
boundary conditions on the inner boundary. Conditions 
at infinity (i.e., fall-off) for are the standard ones, 
that is /3" — > at inflnity. See, for example where 
weighted spaces have been used and ^3 where a com- 
pactiflcation of the exterior region has been employed. 
With this fall-off condition, the kernel is always trivial 
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because there are no conformal Killing vectors which de- 
cays to zero at infinity. However, in the presence of a 
conformal Killing vector normal to the boundary, equa- 
tion H4.21|l still plays a role. In this case it relates some 
asymptotic components of the solution with a boundary 
integral; this is the analog of the restrictions studied in 
|l7ll36l l3^ for the momentum constraint. 



out controlling the sign of this function. We have shown 
that this can be achieved imposing a Neumann (oblique) 
boundary condition on the radial part of /3" and Dirich- 
let conditions on its tangential parts. Using the theory 
of elliptic systems, we have shown that these boundary 
conditions are well posed. 



V. CONCLUSIONS 

In this paper, we have explained why a Dirichlet con- 
dition on /3° may be potentially problematic for solving 
the Hamiltonian constraint. This is essentially because 
with a Dirichlet condition it is not possible to control the 
sign of Kabs'^s^ at the boundary, since this function de- 
pends on normal derivatives of /3° . To control the sign of 
this function is important for two reason: the first one is 
that for physically interesting solutions, i.e. ones which 
contain marginally future trapped surfaces, on maximal 
slices KabS°'s^ is always non-positive. The second one, 
is that with a definite sign of KabS°'s^ it is possible to 
prescribe a priori conditions that guarantee that the so- 
lutions of the non linear equations will always exist. It 
seems to be very hard to obtain such conditions with- 
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